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Abstract 

We present an implementation of the Linked Dipole Chain model for deeply 
inelastic ep scattering into the framework of the Ariadne event generator. 
Using this implementation we obtain results both for the inclusive structure 
function as well as for exclusive properties of the hadronic final state. 



*The original publication was based on results from an implementation containing an error. In 
this revised version this error has been corrected, some of the beyond leading-log assumptions have 
been revised and so have some of the results. 



1 Introduction 



With the HERA colhder, a new kinematical regime has been opened up for studying 
deeply inelastic ep scattering (DIS) on partons carrying a very small momentum 
fraction x of the proton. Much theoretical and experimental effort has been made 
to increase our understanding of the dynamics of such small-x partons, and much 
progress has been made, although many problems still need to be solved. On the 
theoretical side, a major issue is how to handle the resummation of large logarithms 
of X and in a consistent way, while a big obstacle for the experimental analysis 
has been the lack of theoretically well founded event generators which are able to 
describe the measurements made in this kinematical region. 

The conventional way of describing the QCD evolution of the partons within a 
hadron, is to look at ladder diagrams where an incoming parton with large momen- 
tum and low virtuality undergo successive splittings, thus reducing its momentum 
and increasing its virtuality. The contributions from splittings where the parton 
retains only a small fraction z <^ 1 of its momentum or where its virtuality is in- 
creased by a large factor u, are enhanced by large logarithms of 1/z or z/ respectively, 
therefore such ladders need to be resummed to all orders. 

The so-called DGLAP evolution equations handles the resummation of logarithms 
of the virtuality by summing all ladder diagrams where the virtualities, or the trans- 
verse momenta k_i_, are strongly ordered along the chain, while the BFKL 0] equa- 
tions perform a ln(l/2;) resummation, summing ladders with strongly ordered mo- 
mentum fractions, but unordered in k±. Looking only at the inclusive structure 
function F2, both of these approaches are able to explain the steep rise with 1/x 
measured at HERA [^. This is true also for the very small-x region where the 
DGLAP equations are thought to be unreliable. However, the prediction from these 
approximations rely heavily on the assumption of the input parton distributions 
from which the evolution is started, and the absence of a small-x enhancement in 
the evolution can be compensated with a steeply rising input distribution. 

From the F2 measurement alone it is therefore difficult to estimate the relative im- 
portance of the different resummation approaches, and several suggestions have been 
made to instead look at different details of the hadronic final state to get a better 
understanding of the dynamics of QCD evolution Neither the BFKL or DGLAP 
equations are, however, suitable for describing non-inclusive event properties since 
that may destroy cancellations between real and virtual diagrams which are essential 
for the different approaches. In contrast, the so-called CCFM evolution equations 
are designed to explicitly describe exclusive final-state properties by very carefully 
handling interferences between initial- and final-state splittings in the ladder, based 
on an angular ordered description. It can be shown to reproduce both the BFKL 
and DGLAP equations in their respective regions of validity. 
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To properly analyze measured properties of the hadronic final states found at HERA, 
it is important to have event generators which reproduce these properties to a satis- 
factory level. Unfortunately this has not been the case. The conventional generators, 
such as HERWIG and Lepto [0, are based on leading-log initial state parton 
showers derived from the DGLAP equations, and predict much too small partonic 
activity in the direction of the incoming proton (hereafter also referred to as the for- 
ward direction) for small-x events This is expected, as the cascades are strongly 
ordered in transverse momenta, and are therefore limited by the smallness of in 
these events. In contrast, the dipole cascade implemented in the Ariadne [|10| pro- 



gram, where the generated partons are unordered in transverse momenta, is able 
to describe the final-state properties in the proton direction quite well. From this 
one may suspect that the resummation in the BFKL equations indeed are impor- 
tant, but since there is no clear relationship between the semi-classical soft radiation 
model |TT| in Ariadne and the BFKL equation besides the fc_L non-ordering, no firm 



statement can be made. 

It is therefore important to construct an event generator implementing eg. the 
CCFM evolution, which can make reliable prediction about exclusive properties of 



the hadronic final state. Attempts has been made in this direction ||T2[, but several 
difficulties have been encountered. A major obstacle with the CCFM equation is the 
presence of the so-called non-eikonal form factor which makes any implementation 
extremely inefficient. 

Recently a reformulation of the CCFM equation has been proposed [|13|. In this. 



the Linked Dipole Chain (LDC) model, the division between initial- and final-state 
radiation diagrams is redefined using the colour dipole cascade model (CDM). After 
this redefinition, the non-eikonal form factor drops out, which allows for a more 
simple implementation in an event generator. 

A first attempt to construct such an event generator is presented in this paper. 
Although the LDC model is well suited for implementation in a Monte Carlo program 
there are a number of problems to be resolved. One problem is that both the CCFM 
and LDC models only deal with purely gluonic ladders while, to make a complete 
event generator, all types of ladders should be included. It is also important to 
handle energy-momentum conservation in a sensible way. Another important issue 
is the dependence on the input parton densities and the Sudakov form factors needed 
to regularize the poles in the splitting functions and to conserve the total momentum 
in the evolved parton density functions. 

In section ^ we first recall the main ideas of the LDC model, then in section ^ 
we present the different issues involved in implementing the model in a Monte Carlo 
program. To obtain predictions for the hadronic final state we must first, as described 
in section ^, obtain input parton densities which together with the LDC evolution 
will give a satisfactory description of the inclusive cross section. Some results for the 
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Figure 1: Lepton proton scattering with n perturbative ISB emissions. The emitted 
ISB partons are denoted {qi} and the propagators are denoted {fcj}. 

hadronic final state are then presented in section |^. Finally, our conclusions can be 
found in section |. 



2 The Linked Dipole Chain model 

It is well known that the cross section for DIS events is not describable only by the 
lowest order perturbative terms. Still, it is not necessary to consider all possible 
emissions. A large set of them can be summed over and do, in principle, not affect 



the cross section. In the LDC model 131, the emissions that are considered to 



contribute to the cross section are regarded as Initial State Bremsstrahlung (ISB). 
The description of these emissions is based on the CCFM model which is a leading- 
log approximation of the structure function evolution in DIS. The CCFM model has 
been modified by redefining which emissions should be counted as ISB, resulting in 
a much simplified description. 

To describe final state properties of DIS events, one must also consider Final State 
Bremsstrahlung (FSB). In the LDC model this is done within the framework of the 
Colour Dipole cascade Model (CDM) which has previously proved to give a good 
description of parton cascades in hadronic e'^e~ events and DIS. The general picture 
is that the initial parton ladder builds a chain of linked colour dipoles and that the 
FSB is radiated from these dipoles. 

Let {qi} denote the momenta of emitted partons, {fcj} denote the momenta of the 
propagators (see fig. [^) and z^i be the positive light-cone momentum fraction^ of ki in 
each emission: fc+j = 2;_|_jA;_|_(j_i). According to the CCFM model, the emissions that 
contribute to the cross section (ISB) are ordered in rapidity and energy. Furthermore, 

^Note that throughout this paper we work in the j*p centre of mass system with the proton 
along the positive z-axis, except in section ^ when comparing with experimental data, where the 
7* is along the positive z-axis. 
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there is a restriction on the transverse momenta of the continuing propagator in each 
emission: 

kli > z+iql,. (1) 

The weight distribution, dw, of the initial chains factorizes, with the factors dwi 
given by the following expression (a = Sag/n): 



dw = dwidw2 ■ ■ ■ dwn, 

dwi = — n-^^neiz+i,k^i,q±i) . (2) 
A„e is the so called non-eikonal form factor, given by the expression: 



A„e {z, k±,q±) = exp 



(3) 



In the LDC model, the definition of the ISB is more restricted. Consequently, more 
emissions are summed over and the expression for the weight distribution, dwi, is 
changed. The new restriction is that in each emission, the transverse momentum 
{q±i) of the emitted parton must be larger than the lower one of the transverse 
momenta of the surrounding propagators 

q^i > min (k^i, . (4) 

In fig. ^ we show an example of an emission which belongs to the ISB according to 
the CCFM model but violates the restriction in eq. (^) and is regarded as FSB in 
the LDC model. 

Summing over these emissions, the weight distribution of each allowed emission now 
becomes 

_ dz^i d q±i 

dWi = a 7^. (5) 

z+i Tiqii 

This simplification of the expression for dwi is due to the fact that one can interpret 
the non-eikonal form factor, Ane, as a Sudakov form factor, that is, it is equal to the 
probability of not violating the restriction in eq. @. 

By changing variables to the propagator momenta and integrating dwi over the az- 
imuthal angle (in the transverse plane) it can be written approximately as: 
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Figure 2: CCFM LDC: The emission q is regarded as ISB in the CCFM model but 
in LDC model it is regarded as FSB. This is illustrated (above) with fan diagrams 
and (below) in a log p^— rapidity diagram. Rapidity here is defined as \og{p^/pj_) 
and the dashed lines indicate equal p+ . 

It is instructive to look at three different possibilities for two subsequent emissions, 
specified by different orderings of the transverse momenta of the propagators, num- 
bered fci, k2 and (with the same order as in fig. |l|): 



k±i < k±2 < k±3: This is included in DGLAP and gives the same weight as 
there: dw2dws oc 73- -p-. 

k±i < k±2 > k^^. This gives the weight dw2dw^ oc and resembles a hard 



sub-collision with a momentum transfer i ~ A; 



2 

±2- 



/c_Li > k±2 < k±^: Here dw2dw3 oc 
infrared safe! 



■p--p— , it has no factor tJ— and is thus 



The LDC model, without corrections to the leading log approximation, has previ- 
ously been studied and some qualitative results for the structure functions and final 
state properties in DIS have been presented in refs. [13,15-1^. It is found that 
the LDC model, just as the CCFM model, interpolates smoothly between DGLAP 
and BFKL. The emissions along the rapidity axis can be separated in two phases: 
For rapidities closest to the proton (forward) direction, the ISB chain performs a 
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Figure 3: The decrease of log of the parton density with lnp\ at a certain rapidity, 
for different x-values, for running (solid line) and constant (dashed line) coupling. 

BFKL like motion with a constant mean i^^-flow. At a certain distance from the 
photon end, the transverse momentum of the emissions begins to rise to the photon 
virtuahty, as expected by DGLAP. For the structure functions, a DGLAP behaviour 
(exp(constv^— logx)) is shown for moderate values of x, but for small x-values it has 
the BFKL behaviour. 

A prediction from the BFKL model is that logp^ of the emitted partons along the 
chain would be described by a Gaussian distribution with a growing width oc log 1/x. 
From the result of the LDC model, one can clearly see that this BFKL behaviour is 
indeed present for a constant coupling but not for a running coupling. This is illus- 
trated in fig. ^ where the parton density at a certain rapidity is plotted as a function 
of logp_L for different values of x-Bjorken, for constant and running coupling. The 
Gaussian behaviour is observed for constant coupling, while for a running coupling 
it appears to decay exponentially. It seems though that the exponential decay is 
significant only for events with very small x-values and will probably not be visible 
in currently available data. 



The dipole model |]14| was originally developed for final-state parton cascades from 
a quark-anti quark system. The phase space for gluon emission from a. q — q pair 
is approximately given by the triangular area in fig. The gluons are assumed to 
be radiated from a g — g colour dipole and after each emission, the dipole is split 
into smaller dipoles (fig. Qd), which continue to radiate independently under a p±_- 
ordering condition (shaded area). Also g —>■ q + q splittings have been included in 
this model. The size of the dipole triangle is determined by the total q — q invariant 
mass W. 
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Figure 4: (a) The phase space of parton emission from a q — q pair, (h) After 
the first emission, the phase space is split into two triangle which radiate gluons 
independently, hut under a p±-ordering condition, (c) An event with four emissions. 

The momenta {gj} of the emitted ISB partons in DIS are plotted in fig. Due to 
the ordering in positive and negative hght cone momenta, one can insert the ISB into 
a dipole triangle with a size determined by the photon negative light cone momentum 
(left edge) and the positive light cone momentum of the incoming (non-perturbative) 
gluon (right edge). After doing this, the FSB partons can be emitted in a similar 
way as for the q — q parton shower (fig. Pd). The phase space of the FSB is the 
shaded area in fig. ^a. 

3 The Monte Carlo implementation 

The LDC model has been implemented in a Monte Carlo program and some results 
within the leading log approximation have already been presented. Here we present 
a more complete implementation taking into account some non-leading corrections 
and generating complete events all the way down to the final state hadron level. 
What follows is a step-by-step description of the procedure. 

The basic formula for the evolution of the parton densities is 

xf.{x,Q^) = Tf — [ G,,(x,Q2,a;o,fcio)+ (7) 

Sj{Q^, klQ)5ij{\nx - lnxo)]a;o/oi(a;o, k^^), 

Where Gij{x, Q^,Xo, k^o) is the sum of the weights of all chains starting with a parton 
j at some low scale kj_Q carrying a momentum fraction xq, and ending up with a 
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log(W^ 




Figure 5: (a) ISB emissions plotted in a dipole triangle, (b) The solid circles are ISB 
emissions and the dashed circles are FSB emissions. 



parton i carrying a momentum fraction x being hit by a photon with virtuahty Q^. 
The delta function corresponds to the case of no emissions. G is positive definite, 
and to conserve the total momentum, xfj{Q'^,x)dx = 1, the delta function is 

multiplied by a Sudakov form factor, Sj{Q'^,k'j_Q), representing the probability that 
the parton j with with momentum fraction xq at fc^Q has not split, and thus reduced 
its momentum, when probed at a higher scale Q^. This form factor will be discussed 
in detail below. 



The analytic approximate upper limiting function for G |T^, given by 



G,,{x,Q^xo,kl,)<G{Qykl„x/xo) = ^,h{2Vab) (8) 
a = y/^ilnQ'^/k^Q + liaxo/x), b = \/^\iaxo/x 



provides us with the starting point, and for each generated chain c a number of 
multiplicative weights are calculated Uc = Hiu^c^ so that the correct form of Gij is 
obtained as 

Gij{x,Q'^,Xo,kl^) = uJij{x,Q'^,Xo,klQ)G{Q^/klQ,xo/x), (9) 
with the average weight 



^ = ^E^c (10) 



c=l 

This is how it is done: 
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1. First the x, and the flavour i of the struck quark is chosen using evolved 
parton densities and the standard Born-level electro-weak matrix elements. 
This is currently done within the Lepto program 

2. Then the Xq and flavour of the incoming parton is chosen according to eqs. 
and with k±o as a given parameter. 

3. The number of emissions is chosen from 



-,nl,n— 1 



(11) 

4. The positive and negative light-cone momentum fractions, Zj^ and which 
enters in each emission j is generated according to the ordered integral 

I /a"n,^^5(lnxo + $^ln^,+ -lna:) (12) 

n\[n — l)\ J Zj+ Zj_ ^ 

5. At this point we need to choose the flavours of each link. To do this we 
introduce the standard Altarelli-Parisi splitting functions and preliminary we 
use their approximations, Pi^j{z), in the z —>■ limit 

1 + ^2 

Pq^q{z) = Cp- ~ Cp, 

1 — Z 

z[l — z) z 
Pg^,iz) = Tniz' + il-zf)^Tn, (13) 

z z 

We get the first weight factor as the approximated splitting functions summed 
over all possible flavour combinations. 



a;(°) = 5^n5=^. (14) 



6. We can then use eq. (0) to generate a specific flavour combination according 
to their individual weights. 

7. Next, we generate the azimuthal angles of each emission and construct the exact 
kinematics. The delta function in eq. (0), which handles the conservation of 
positive light-cone momenta, does not take into account the transverse degrees 
of freedom. In particular it would give zero positive light-cone momentum 
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for the struck quark, in the final-state. We therefore modify this delta 

function to exactly conserve the total energy and momentum, effectively setting 
Zn+ by hand to the value needed. However, for some values of the azimuth 
angles this is not possible and we get a weight factor corresponding to the 
allowed integration area A(^j: 

= U~ [ dcP,. (15) 

8. Then we implement the condition that the transverse momenta (which is gen- 
eralized to the transverse mass m± for massive partons) of an emitted parton 
must be larger than the smallest virtuality Vj of the connecting links j and 
j — 1, and that all virtualities must be above fc^Q, giving us the second weight 
factor 

C^(^) = IljQ{m]_j - Vjmm)Q{Vjmin " k]_o) ■ (16) 

We note that the weights are finite even if one of Vj and Vj^i goes to zero, 
and one could imagine replacing the second theta function in eq. (|T6|) with 
©("i^jmax — ^io)- This would reduce the dependency on k±Q and would allow for 
more unordered chains as discussed below. 

9. Now we introduce the running of a^, giving a fourth weight 

^^'^ = n,- , I (17) 

^ln(m^^./A2) ^ ^ 

10. Having obtained the virtualities of the links, we can now correct the splitting 
functions in eq. (0). We get the following cases: 

• fj+i > Vj > Vj-i. Going upwards from the proton side we use P/^._^^/^.(2;j+) 

• Vj^i < Vj < fj-i: Going upwards from the photon side we use P/^^-^^/^ (2;j_) 

• fj+i < Vj > Vj^i. Corresponds to a Rutherford scattering and Zj^ ~ 
Zj- = z. Here we use 2-^2 matrix elements taking into account colour 
connections as explained in Appendix A. 

• fj+i > Vj < Vj-i. Here we use the z ^ limit, Pi-,j{z) of the splitting 
functions, where z = z+ if Vj+i > Vj-i or else z = Z-. 

Note that the colour factor is independent of the ordering of the virtualities. 
Therefor we have to correct for only the kinematical part of the splitting func- 
tion by using reduced splitting functions /'where the colour factor is divided 
out. Here we also introduce the Sudakov form factor, to be discussed below 
and we can write the fourth weight factor 



(4) _ n 



Uj ' ^' - / ' - - ^^-^^^-^^-+1 - ^ (18) 
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where Sj{vj_i, Vj, f j+i) is Sf._-^{vj,Vj^i) or Sf.^-^{vj, fj+i) depending on whether 
the virtuahty is going up or down. 

Pq^q and Pg^g both have poles as z — > 1, corresponding to emission of low- 
energy gluons. Typically these should be counted as final-state emissions, but 
to be sure to avoid divergences we introduce a cutoff z^^^ = 0.5. See also the 
discussion of double counting below. 

Note also that we use mass less splitting functions, and the production of heavy 
quarks is only suppressed by the phase space. This should be improved in the 
future. 

11. The final weight factor is introduced to correct the emission closest to the 
photon, in the cases where f„ > fn-i, to reproduce the exact 0{aas) matrix 
element as given eg. in ref. ||18||: 



(5) _ ■M.jQ'^jX, Zn+ , Zn- ) , s 

12. The generated chain is now kept with a probability ujc = '^^i^c \ where W 
is a scale factor to avoid probabilities larger than one. There is in principle 
nothing preventing weight larger than one, but they turn out to be very rare. 
Nevertheless, it may happen, and it is important to check that W is large 
enough so that the results are not infiuenced by this. Chains with u > 1 may 
optionally be saved and retrieved again when an event with the same fiavour 
and similar x and is requested. A chain will then be used on the average uJc 
times, each time with different final-state cascade and hadronization. Below 
we have used W = 1 giving less than 0.1% events with weight larger than one. 

13. To prepare for the final-state dipole radiation the emitted partons must be 
connected together and form dipoles. In the case of quark links, this is straight 
forward, the incoming quark is simply connected to the first emitted gluon and 
so on until the struck quark. In the case of gluon links, there are two colour 
lines, and a radiated gluon can belong to either of these. This choice is done 
completely at random. The connection between the colour line of the incoming 
parton and the proton remnant is handled in the same way as in the default 



soft radiation model of Ariadne |[T0|] . 



One could imagine using other methods for determining the colour-fiow. One 
suggestion is to use the colour-fiow which minimize the total string length]^. 
One could also consider colour-reconnections, following eg. the model already 
implemented in Ariadne , possibly giving rise to large rapidity gaps among 
the final-state hadrons. 



^As defined eg. by the A measure of ref. 
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14. The constructed dipoles can then radiate more final-state gluons in the phase- 
space hmited by the virtuahty and the positive and negative hght-cone mo- 
menta of the hnks in the chain as in fig. Note that also in the no-emission 
case corresponding to the delta function in eq. some radiation is allowed 
within the triangular area defined by g+ < —Q+ and g_ < Q-. In the Breit 
frame, this is just the area of allowed FSB in a, e^e' qq event with centre 
of mass energy Q. 

15. Finally the final state dipole chains are hadronized according to the Lund string 



fragmentation model as implemented in Jetset [pl|] . 

This concludes the description of the actual implementation. But before we can start 
producing events we have to fix the parameters involved. These are A, k±o, W, z^^t 
and the input parton densities xofoj{xo, fc^o)- ^ is not really a physical parameter, 
and should be set large enough so that the result no longer depend on it. It would be 
natural to take A and k±o to be the values which have been tuned for the final state 
dipole cascade to reproduce LEP data, Alep and /cxolep- One could, of course, use 
an increased cutoff k±o > ^xolef in the ISB. This would mean that more emissions 
would be moved from the initial to the final state, which would continue emitting 
down to A;_lolep- We therefore expect the final result to be fairly stable w.r.t. such 
variations as long as /c_lo is not too large. 

There is an additional complication with a large k±o for the cases where the virtuality 
drops below k±Q somewhere along the chain, causing a zero weight in eq. (0). For 
the total cross section, this does not matter, as such fiuctuations are included in the 
input parton densities at a lower xq, corresponding to the momentum fraction of the 
link closest to the photon which is below the cutoff. For the final state, however, it 
means that we are excluding some radiation close to the direction of the incoming 
hadron. 

As discussed above, one could replace the second theta function in eq. (|16|) with 
©('fjmax — ^io) since the then would still be in the perturbative region. In this 
way the result would be less sensitive to variations of k±Q, as the perturbative system 
on the proton side of the sub-cutoff link would still be generated. This would not, of 
course, solve the problem altogether as one can imagine chains where two or more 
consecutive links are below the cutoff. 

In step |l^ we have replaced the Nc/z pole, which is used in the emissions of the 
original leading log LDC model, with the standard Altarelli-Parisi splitting functions. 
This should be a sensible way of including some sub-leading effects, as long as the 
splitting functions are regularized in a correct way. In each emission, two particles 
are produced and one vanishes. For the parton distributions, this means we must 
subtract and add partons correspondingly. 

A simple way of treating this double counting problem is to add only one of the 
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produced partons assuming that the mother parton is not affected by the emission. 
The choice of which of the two partons to add is not trivial. One way is to introduce 
a cut-off Zcut = 0.5 allowing only z < Zcut- This is a good choice in the g ^ g 
splitting since the gluon with z > 0.5 is more similar to the mother gluon. For the 
other emissions, which involve both quarks and gluons, it is more important to make 
the choice which leads to a better approximation of the quark distributions. This 
can be done by allowing all q g splittings, forbidding all q —>■ q splittings and 
allowing g ^ q splittings only if the quark (or anti-quark) interacts directly with the 
photon. 

In a more sophisticated treatment, the subtraction of a parton is done with Sudakov 
form factors corresponding to the probability for the partons not to vanish before the 
splitting can occur. In this way, we can take into account that e.g. the possibility for a 
gluon to spht into a quark anti-quark pair with low virtuality reduces its contribution 
to emissions with higher virtualities. The suppression factor becomes an exponential 
of an integral over the splitting 



l^'p,^,{z)dz I 



^, = exp -/ P,^,iz)dz I ^^-^ 



(20) 



The region of integration corresponds to the region of allowed emissions. Here we 
only use an approximate form. A more thorough investigation will be presented in 
a future paper. 

The lower limit on is given by the lowest virtuality of an emission step multiplied 
by a fudge factor, e^, to account for the suppression of the emission probability for 
small g^-values (eq. (|)). The value of 5 = 0.4 that is used has proven to be an 



effective cut-off in the leading log treatment of the LDC model |T^,|T3- The upper 
limit is set to the highest virtuality of an emission step. One could imagine having 
a higher, or a z dependent limit, but for simplicity we only use the virtuality in this 
publication. 

The Sudakov factor due to the g g splitting depends on the choice of Zcut- For 
Zcut = 0.5, there is no double counting to correct for and for larger Zcut the integration 
of the splitting function is in the region 0.5 < z < z^ut- H q ^ q splittings are allowed, 
quarks are suppressed with a Sudakov factor where the integration region is given 
by < z < Zcut- 

The situation is quite different when we are interested in the final state properties. 
Here, the Sudakov form factors are not as important since they for most events 
roughly give an overall factor YYj=i S{vj-i,Vj,Vj+i) ~ S{k'j_Q,Q^) which does not 
have an influence on the relative contributions of different final states. On the other 
hand, disallowing some of the initial state emissions to reduce double counting has a 
large effect on the final state properties since it reduces high p± emissions. This is a 
problem except for gluon emissions with z > 0.5 in regions with ordered virtuality, 
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Figure 6: In the case where the virtuality is increasing, a splitting with z larger than 
0.5 corresponds to an emission of a gluon q which is inside the phase space where 
subsequent FSB emissions is allowed. Such splittings should therefore he removed 
from the ISB chains to avoid double counting. 



since these gluons can be treated as final state emissions as shown in fig. Conse- 
quently, we would expect it to be a good approximation for the final state generation 
to skip the Sudakov form factors and to allow some of the sphttings which lead to 
double counting. 

It is clear that the result is very dependent on the non-perturbative input parton 
densities, which are basically unknown. If eg. the gluon density is very divergent at 
small X, the Xq chosen from eq. will tend to be small, limiting the total phase 
space available for radiation Ay oc Iiixq/x. The input parton densities can, however, 
be constrained somewhat from the total cross section, and we can parametrize them 
and make a fit of the parameters to eg. F2 data at different x and . 



4 Fitting the input parton densities 

From eq. we can write the leading order expression for F2 as 

F2 = Y.^]Y.f — ^ G,,(x,Q^xo,fcio)+ (21) 

~7n A Jx ^0 



SjiQ'^, klQ)6ij{\nx - Inxo) ] Xofoj{xo, kl^). 



For given values of x, Q^, xq, k^^Q, i and j, we can calculate Gij{x, Q^, xo, fc^g) from 
steps ^ through 11 in the previous section using eq. (P). For any given parametriza- 



tion of the input densities it is then possible to calculate F2 and compare with 
experimental data. 

In principle one could also fit to other data, such as prompt photon and jet production 
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in hadron-hadron collisions. This is not possible in our current implementation which 
only gives the evolved densities for quarks and anti-quarks. This means that the input 
gluon distribution is only constrained indirectly from the dependence of F2, and 
in this paper we only make a very crude fit using only four different parameters for 

The input densities are parametrized as 

xh,[x) = A.x'^^iX-xfK (22) 

For the valence distributions u^ix) and d^{x) we use the same form, with = 3 
leaving free and using the normalization 

/ Uy{x)dx = 2, / dy{x)dx = 1 (23) 
Jo Jo 

to fix Auv and A^v- For the gluon distribution, Pg = 4: while ag and Ag are left free. 
All the sea-quarks distributions have the same form with = 4, leaving as free and 
setting 2As = 2As = Au = A^ = Ad = A^ so that the total momentum 

f dxy^xh,{x) = l (24) 
Jo 

is conserved. 



We only use data from proton F2 measurements from HI ||2^, ZEUS [^], NMC ||24 



and E665 p5[ without allowing for any normalization uncertainty factors. We use 
only data for > 1.5 GeV and x < 0.5 to ensure a reasonable length of the 
evolution. We then make six sets of fits using different options in the generation of 
the G function: 



A LDC default: fc^o = 0.6 GeV, A = 0.22 GeV. 

B DGLAP: As for A but only allow chains with monotonically increasing virtu- 
alities of the links from the proton side. 

C DGLAP': As for B, but chains where the virtuality of the link closest to the 
virtual photon is larger than are permitted. We use this as a kind of higher- 
order corrected DGLAP evolution although, of course, not equivalent to NLO 
evolution. 

D As for A but k±Q = 1 GeV, to check the sensitivity to this cutoff. 

E As for A but without the Sudakov form factor. Instead Pq^q{z) is set to zero 
and Pg^q{z) is nonzero only in the splitting closest to the photon. 
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Figure 7: The fitted F2 as a function of x for different strategies and for two different 
values of Q"^ , compared with data from 

F As for A but (3g = (3s = to check the sensitivity to the fit parameters. 

G As for A but only fitting to F2 data with x < 0.1, to reduce the sensitivity to 
the step size 5\n.xo/x = 0.2 used when integrating eq. (PD , and to the high-x 
form of the input parametrization. 

H As for A but allow the virtuality of some links to be below kj_Q as long as the 
largest virtuality of two consecutive links always is above kj_Q. 

The results of the fits are presented in figs. ^ and ^. For the default case, the fit is 
quite acceptable. We note in particular that the fitted input gluon density is slowly 
decreasing with 1/x, although we must keep in mind that the gluon distribution is 
only indirectly constrained. 

For the DGLAP case the fit is much worse. The number of allowed ISB emissions is 
here strongly restricted, especially for small and x. This results in much slower 
evolution which forces the input densities to rise with 1/x. In fit C, where the link 
closest to the virtual photon is allowed to be above Q^, the fit is on the other hand 
again quite acceptable. The input gluon distribution is no longer rising with 1/x, 
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Figure 8: The fitted input gluon density as a function of x for different strategies. 



on the other hand it is as strongly decreasing as in fit A. It is known that F2 can be 
fitted using conventional DGLAP evolution with a valence-like flat input distribution 



at small input scales, close to the one used here, as in the GRV parametrizations ||26|| . 
Such fits can, however, not be directly compared to this one, as we here have less 
parameters. But we can conclude that the ISB chains with unordered k± do play a 
role in our case, although most of the effect can be obtained allowing only for one 
'stepping down', closest to the photon in the chain. 

Increasing the input scale to 1 GeV in fit D makes a big difference particularly at 
small X as seen in figs. 0d and §3. Also here the number of allowed ISB emissions is 
strongly restricted and again the input gluon is forced to increase strongly with 1/x. 
This should come as no surprise. We expect, however, that when we below study 
the hadronic final states, the result should be less sensitive to the input scale used. 

The importance of the Sudakov form factor is apparent in figs. |^ and ^b, comparing 
fits A and E, especially for the input gluon distribution. Also in figs. and ^ we 
show the fit H, where some propagators below the cutoff is allowed. The reproduction 
of F2 does not change much, but we see that the input gluon decreases slightly faster 
with 1/x than for the default fit A. 
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The fits F and G in figs. |^ and ^ show how sensitive the fit is to the input distri- 
bution is at large x. Changing f3g and f3s from 4 to 5 does not influence the fit very 
much and neither does the omission of the data points at high x, although in both 
cases the input gluon is shifted somewhat to higher x. 



5 Results for the hadronic final states at HERA 

The results for the hadronic final states will depend on the input densities used. If 
eg. the input densities are increasing with 1/xo, chains starting with low xq will be 
favoured and the length of the chains will be shorter resulting in fewer emissions 
and less activity in general due to the reduction of the available phase space, in 
particular close to the direction of the incoming hadron. Also a smaller input gluon 
density will result in fewer chains initiated with a gluon, which has higher charge than 
an incoming quark, and again the probability of emissions will become smaller. In 
addition, the non-perturbative hadronization is smaller if only one string is stretched 
between the perturbative system and the hadron remnant in the case of an incoming 
quark. 



To study the hadronic final states we use the HZTOOL package developed jointly 



by HI, ZEUS and theoreticians for comparison between event generators and pub- 
lished experimental data. We have selected six different distributions which have 
been shown to be sensitive to details in the models used in Monte Carlo event gener- 
ators. The distributions presented in fig. ^, which are all measured in the hadronic 
centre of mass system, are as follows. 

(a) The i^^_L-fiow as a function of the pseudo rapidity for two bins in x and Q^, one 
with low X and low Q^, and one with moderate x and in (b) The large 
amount of E± in the forward direction was previously claimed to be a good 
signal of k± non-ordering in the ISB, but it has been shown that this effect can 
also be obtained by the introduction of additional non-perturbative effects. 

(c) The so-called seagull plot with the average fc^ as a function of Feynman-x 
[ p^ , which at eg. EMC ||2^ was shown to be difficult to reproduce with event 
generators. 

(d) The ^^-distribution of charged particles in a forward pseudo rapidity bin. This 



was recently proposed [§, 0] as a new signal for perturbative activity in the 
forward region indicating k± non-ordering: a h.igh.-k± tail would be difficult to 
reproduce by non-perturbative models, where such tails would be exponentially 
suppressed. 

In was also shown that the pseudo rapidity distribution of charged particles 
with k± > 1 GeV also could be a good signal for k± non-ordering. 
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Figure 9: Comparison of the default LDC A model with HERA data as given in [W^ 
and other event generators. The distributions are as follows: (a) The transverse 
energy flow as a function of pseudo rapidity for events with 0.0001 < a; < 0.0002, 



5 < qyCeV < 10 



(b) As (a) but for events with 0.003 < x < 0.01, 20 < 



Q"^/ CeV^ < 50. (c) The average squared transverse momentum of charged particles 
as a function of xf /p^. (d) The transverse momentum distribution of charged 
particles in the pseudo rapidity bin 0.5 < r/ < 1.5 for events with 0.0002 < a; < 0.0005, 
6 < Q"^/ Ge\^ < 10 J^. (e) The pseudo rapidity distribution of charged particles with 
a transverse momentum larger than 1 GeV for the same kinematical bin as in (d). 
(f) The two-jet ratio R2{x) as a function of x ^3^] . All measurements were made 
in the hadronic centre of mass system and only events without a large rapidity gap 
were included. The full line is LDC A, long-dashed is Ariadne 4-08 with default 
parameter settings, dotted is Lepto 6.4 with default parameter settings and short- 



dashed is the same but with SCI and the special sea-quark remnant treatment ^SSj] 
switched off. 
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Figure 10: Comparison between different LDC strategies and data corresponding to 
fig. ^e. In all cases the full line is the default LDC A strategy. In (a) the dashed 
line is LDC B using only DGLAP-like chains and dotted is LDC C, i.e. the same 
but allowing the virtuality of the link closest to the photon to be above . In (b) 
the long-dashed line is LDC D with k±Q = 1 CeV, the short-dashed is LDC with 
increased cutoff, z^-y^f, in the splitting functions and dotted is LDC H, allowing some 
propagators below the cutoff. In (c) the dashed line is LDC E without Sudakov form 
factors and disallowing all q ^ q and most g q splittings in the parton density 
fit and the dotted, LDC Aq, is the same as LDC A but using Sudakov only in the 
parton density fit. Finally in (d) the dashed line is LDC F using a different form of 
the input densities at large x and the dotted is LDC G restricting the fit to F2 data 
with X < 0.1. 
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(f) Result for the two-jet ratio R2 was recently reported and showed large 
differences between the currently used event generators. 

In fig. ^ we see the results from LDC with default settings and using fit A, labeled 
LDC A, compared with data and with the results for Ariadne using the default soft 
radiation model, with Lepto with and without the additional non-perturbative (soft 
colour interactions (SCI), and perturbative-like treatment of remnants in the case of 
sea quarks) assumptions presented in [Q. For the Ariadne and Lepto models we 
confirm previous results where Ariadne, which until now was the only model im- 
plementing fcx non-ordering, reproduces the data very well while the DGLAP-based 
Lepto has difficulties, especially without the additional non-perturbative models. 
We find that the result for LDC is quite acceptable, although not reproducing data 
as well as Ariadne. 

To compare different LDC strategies, we show in fig. [1^ only the number of charged 
particles with transverse momentum larger than 1 GeV as a function of pseudo 
rapidity for small x and (fig. [l^). The effects on the other distributions in fig. | 
are very similar. 



In fig. p!0| a we see the results for LDC when restricting to DGLAP-like chains. In this 
case, corresponding to the lines marked LDC B, the result is very poor as expected. 
Allowing the virtuality of the link closest to the photon to be above as for LDC C, 
makes things much better and only slightly worse than the default LDC A. Naively 
one may expect this to give the same result as Lepto which uses the exact 0{as) 
Matrix Element for the emission closest to the photon, also allowing the first link to 
have a virtuality larger than Q^, and which adds on parton-showers a la DGLAP on 
such configurations. But in LDC, even though no ISB emissions are allowed between 
the highest virtuality link and the photon, there is still a resummation of diagrams 
which are then replaced by FSB emissions. Lepto, however, uses the 'bare' matrix 
element and does not include any resummation. 

For the line marked LDC D in fig. pi]|b, the cutoff in k± for the ISB is set to 1 
GeV. In the previous section we saw that the parton density functions in this case 
became dramatically different. For the final state, however, the reduction of ISB is 
compensated by final state dipole emissions, which are allowed in the whole rapidity 
range below /c_l = 1 GeV down to the cutoff fitted to LEP data, /c^olep = 0.6. The 
effect of increasing z^nt is also shown in fig. p!0| b for the line LDC A^. The dependence 
on this cutoff is small, which is expected as most of the emissions with 2; > 0.5 are 
counted as FSB as explained in fig. ^ above. Also in fig. pi^ j is shown the effects 
of allowing the virtuality of some links below the cutoff. Again the differences are 
small. 

In fig. |TD|c we see the effect of different regularizations of the splitting functions. The 
line LDC E uses the fit E in the previous section, where aX\ q ^ q and most of the 
g —>■ q splittings are disallowed and all Sudakov form factors set to 1. For the final 
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state generation, all splittings are again included, and we see a clear enhancement 
w.r.t. the default LDC A strategy which includes Sudakov form factors. The line 
LDC Ao uses the same parton densities as the default strategy, but excludes the 
Sudakov form factors when generating the final states. The main effect of the form 
factor is to scale down all weights, and since in the total number of events is fixed, 
the effects on the final state is small. The fact that LDC E is as different is then 
mostly due to the difference in the input parton densities. 

Finally in fig. |TD|d, we show the effects of using different fitting procedures for the 
input parton densities and we see that the differences are small. 

6 Conclusions 

We have here presented the first implementation of the Linked Dipole Chain model 
in an event generator. Being based on the CCFM formalism it represents one of the 
first attempts to correctly describe the details of the hadronic final states in small- a; 
deep inelastic scattering to leading-log accuracy. 

The LDC model was originally formulated only ioi g —>■ g splittings in a strict leading- 
log approximation. Going from this to a full event generator is not trivial, and we 
have here described how we implement massive quarks splittings, sub-leading correc- 
tions, convolution with input parton densities, energy and momentum conservation, 
Sudakov form factors, final-state radiation and hadronization. 

Our implementation still suffers from some uncertainties. One is the input parton 
densities which are poorly constrained because only F2 data can be fitted, where 
the gluon only enters indirectly. Another issue is the uncertainty in how to deal 
with the final state of chains where a link drops below the cutoff, surrounded by one 
perturbative system on each side. But from the results in the previous sections, the 
main uncertainty is the regularization of the splitting functions and the conservation 
of total momentum in the evolved parton density functions. The correct way to treat 
this is with Sudakov form factors. But in the current implementation these are only 
treated in an approximate way and they need to be examined in more detail in the 
future. 

Despite these uncertainties, the result presented here allows for some conclusions. 
Compared with the Lepto event generator, which is based on DGLAP evolution, 
LDC is clearly better in describing the hadronic activity in the forward region at 
HERA at small x. The description is not perfect, however, and it seems that LDC 
is still underestimating the perturbative activity in the forward region. 

We have also shown the relative importance between DGLAP-like chains with mono- 
tonically increasing virtualities and the unordered chains in the full LDC model, and 
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found that unordered chains indeed are very important, but that at HERA, most of 
the activity can be attributed to DGLAP-hke chains where the hnk closest to the 
photon is allowed to be above . 

The main goals for the future developments of the LDC generator is to investigate 
the exact form of the Sudakov form factors, and to include a treatment of hadron- 
hadron collisions, to get a better constraint on the input parton densities. Also 
further studies of the final state of chains with two or more perturbative systems 
connected with sub-cutoff links ought to be done. This is especially interesting in 
connection with the large fraction of rapidity-gap events found at HERA. 

Despite the shortcomings of the current implementation, we feel that the LDC event 
generator may become a very important tool for understanding the small-x hadronic 
final states at HERA. 
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Appendix A: Colour connected matrix element cor- 
rections 

We will here describe the splitting functions that are used for the local sub-collisions 
(see point in section These are derived from the corresponding colour connected 
2 — s> 2 QCD matrix elements presented in [^ . 



In order to use matrix elements for deriving splitting functions, it is necessary to 
take colour connections into account. We can illustrate this statement with the 
matrix element for the g + g q + q process. Suppose we are in a frame where 
the incoming gluons are moving head on and parallel to the longitudinal axis and let 
C+(-) be the positive (negative) light-cone momentum fraction for the quark and C+(-) 
be the same for the anti-quark. Because of symmetry between the quark and the 
anti-quark it is clear that the matrix element must be symmetric in these variables: 
|A^P(C+(_), C+(-)) = l-^P(C+(-)) C+(-))- III program, the colours are connected 
randomly so that e.g. the quark has the same probability to be connected with the 
gluon coming in from either the proton or the photon side. This would mean that 
the two diagrams a') and M) in fig. |ll| would be equally probable. 

This is because the matrix element is calculated with an averaging of the colour 
states of the gluons. The corresponding colour connected matrix element, |A1ccP 
is on the other hand calculated with the assumption that the colour states of the 
incoming gluons are known. This means that we get two separate contributions, one 
from the case that the quark is colour connected with the gluon coming in from one 
side, say the proton side, and one contribution from the case that it is connected 
with the gluon from the photon side. Each of these contributions is non-symmetric 
with respect to the quark and anti-quark variables, while the sum of them still is 
symmetric. 

A similar procedure is used in the 0{aas) matrix element correction (see point ^ in 
section]^), to separate the '-/g matrix element into two contributions corresponding 
to quark and anti-quark scattering respectively as described in ref. ||34|| . 

For a colour connected matrix element \A4cc\'^{s,i,u), the corresponding splitting 
function is given by the formula 

P{z) (X z{l- z) ■ \Mcc\'^{l,z,l- z) with z = z+^z_. (25) 

The splitting functions f^j^i^z) where the colour factor is divided out and with the 
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Figure 11: Feynman vs. colour connection diagrams for the g + g Q + q sub- 
process. If we do not average the colour states, the diagrams a ') and b ') give separate 
contributions. The contribution from diagram b ') should be suppressed. 



successive flavours ijk for the propagators, become 



- (l-zf[z' + {l-zf] 
KJ^) = f9j^) = (26) 



ggq' 

P ( 

gqg^ 



999 



z) = {l-zf[z' + {l-zf] 

1 

Yz 



999(')=UZ) = l[l+(l-.f]^ 



= ^^{^l^£^.eio.5-z)^2'-l[i-z^z^r 



The ^-function in the first term of f^gg is a cut-off to prevent divergences. This is 
needed since the Sudakov form factors that arc used regularize only regions with 
ordered virtualitics. The splitting functions f^^^ and f^^^ are set to zero since they 
correspond to the same {g + q g + q) scattering as the splitting functions f^^^ and 
fggg ^^^'^^ have been chosen to take the whole contribution. 

Here we have summed the contributions from different colour connections for each 
fiavour combination. The choice of colour connections has an effect on the multiplic- 
ity and transverse momentum distributions at the hadronic level. Therefore, we will 
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in the future use this information also to choose the colour connections for the local 
sub-collisions. 
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